Correlation functions for a di-neutron condensate in asymmetric nuclear matter 
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Recent calculations with an effective isospin dependent contact interaction show the possibility 
of the crossover from superfluidity of neutron Cooper pairs in ^5*0 pairing channel to Bose-Einstein 
condensation (BEC) of di-neutron bound states in dilute nuclear matter. The density and spin 
correlation functions are calculated for a di-neutron condensate in asymmetric nuclear matter with 
the aim to find the possible features of the BCS-BEC crossover. It is shown that the zero-momentum 
transfer spin correlation function satisfies the sum rule at zero temperature. In symmetric nuclear 
matter, the density correlation function changes sign at low momentum transfer across the BCS- 
BEC transition and this feature can be considered as a signature of the crossover. At finite isospin 
asymmetry, this criterion gives too large value for the critical asymmetry ~ 0.9, at which the BEC 
state is quenched. Therefore, it can be trusted for the description of the density-driven BCS-BEC 
crossover of neutron pairs only at small isospin asymmetry. This result generalizes the conclusion 
of the study in Phys. Rev. Lett. 95, 090402 (2005), where the change of sign of the density 
correlation function at low momentum transfer in two-component quantum fermionic atomic gas 
with the balanced populations of fermions of different species was considered as an unambiguous 
signature of the BCS-BEC transition. 

PACS numbers: 21.30.Fe, 21.60.-n, 21.65.Cd 
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I. INTRODUCTION 



Few-body correlations play an important role in low- 
density nuclear systems such as expanding nuclear matter 
in heavy ion collisions, the tails of nuclear density distri- 
butions in exotic nuclei, or nuclear matter in the crust of 
a neutron star. Feasible low-density phenomena include 
a-particle condensation [l|, the formation of two-body 
(deuteron) 0] and three-body (triton) 0] bound states, 
and their possible coexistence/competition [1]. Of spe- 
cial interest also is the crossover from superfluidity of nu- 
cleon Cooper pairs to Bose-Einstein condensation (BEC) 
of two-body bound states in dilute nuclear matter. The 
transition from large overlapping Cooper pairs to tightly 
bound pairs of fermions can be described on the basis of 
BCS theory, if the effects of fluctuations are disregarded. 
First, this transition was studied in excitonic semicon- 
ductors 0, ordinaiT superconductors Q and in an at- 
tractive Fermi gas Although the BCS and BEC hm- 
its are physically quite different, the crossover between 
them was found to be smooth within the BCS theory. In 
nuclear matter, the BCS-BEC crossover is usually men- 
tioned in respect to the formation of Bose-Einstein con- 
densate of deutcrons at low densities, which is caused 
by strong neutron-proton (np) pair correlations in the 
'^Si Di pairing channel [3, l8|. During this transi- 
tion the chemical potential changes sign at certain crit- 
ical density (Mott transition), approaching half of the 
deuteron binding energy at ultra low densities. Despite 
the fact that in the BCS regime even small isospin asym- 
metry effectively destroys a np condensate [1, 0, [l3| , at 
low densities Bose-Einstein condensate of deuterons is 
weakly affected by an additional gas of free neutrons even 
at very large asymmetries [ll|, [l2| . 



There is a rather general reason to expect that di- 
neutron correlations in low-density nuclear matter should 
be of significance as well. As is known, the bare nucleon- 
nucleon interaction in the ^Sq channel leads to a virtual 
state around zero energy characterized by the large neg- 
ative scattering length a„„ ~ — 18.5fm. This implies a 
very strong attraction between two neutrons in the spin 
singlet state that may lead to the formation of compact 
bound pairs of neutrons in low-density nuclear systems. 
Recent theoretical studies confirm this suggestion. In 
particular, the spatial structure of the neutron Cooper 
pair in low-density nuclear and neutron matter was an- 
alyzed with the use of the bare GSRS force and the ef- 
fective Dl Gogny interaction [3]. It was shown that the 
size of the Cooper pair becomes smaller than the average 
inter-neutron distance at a wide range of neutron densi- 
ties Qu/qo ~ 10~* — 10~^, independently of the type of 
the interaction {go = 0.16 fm"'^ being the nuclear satu- 
ration density). The neutron pair wave function demon- 
strates the BCS-BEC crossover behavior, being spatially 
extended and of the oscillatory type at moderate den- 
sities and possessing the well localized single peak at 
low densities. The spatial structure of the wave func- 
tion of two neutrons was studied also in the halo nucleus 
^^Li [l3|. The behavior of the neutron pair at different 
densities was simulated by calculating the pair wave func- 
tion at several distances between the core nucleus and the 
center of mass of the two valence neutrons. Analogously 
to the case of infinite nuclear matter, the same BCS- and 
BEC- like structures in the pair wave function were found, 
corresponding to the normal and low density regions, re- 
spectively. In Ref. p^ . the di-neutron correlations were 
studied on the base of a new effective density dependent 
contact pairing interaction. An important feature is that 
this contact interaction contains the isospin dependence 



2 



and was adjusted to reproduce the pairing gap in sym- 
metric nuclear matter and neutron matter, obtained in 
the microscopic G-matrix calculations with the realistic 
Argonne vis nucleon-nucleon (NN) potential [l^. This 
study shows that for the screened pairing interaction a di- 
neutron BEG state is formed in symmetric nuclear mat- 
ter at neutron densities gn/ga ^ 10~^ and there is no 
evidence to the appearance of the BEG state in neutron 
matter. 

Since a two-neutron system in free space has a vir- 
tual state around zero energy, then the medium polar- 
ization effects should play the major role in the possi- 
ble formation of the di-neutron BEG state in low-density 
nuclear systems. However, the influence of the nucleon 
medium on the di-neutron correlations seems to be dif- 
ferent in neutron matter and in symmetric nuclear mat- 
ter. Most of the recent quantum ma ny-body calculations, 
e.g., based on the Brueckner theory [ia[l3i O'" renormal- 
ization group approach [isj . show the reduction of the 
pairing gap in neutron matter by a factor of 2-3 due to 
the screening effect of the pairing interaction. The only 
exception is the variational quantum Monte Garlo cal- 
culations of Ref. [l^, showing a small influence of the 
medium polarization on the neutron ^5*0 pairing. At 
the same time, calculations based on the Brueckner the- 
ory show that in symmetric nuclear matter the medium 
polarization of the interaction enhances the di-neutron 
pairing correlations due to the proton particle-hole exci- 
tations p^ . The same evidence in favor of the medium 
enhancement of the neutron ^ So pairing was obtained for 
finite nuclei, where the pairing interaction becomes more 
attractive because of the surface vibrations [13, [2l|, . 
Thus, one can conclude that the conditions for the for- 
mation of a di-neutron BEG state are more favorable in 
a nuclear environment rather than in neutron matter. 

The recent upsurge of interest to the BEG-BGS 
crossover is caused by finding the BGS pairing in ultra- 
cold quantum atomic gases [23l . [13 ■ Ultracold atomic 
gases provide an experimental play-ground for testing 
pairing phenomena due to the possibility to control the 
inter-atomic interactions via a magnetically-tuned Fes- 
hbach resonance. In fact, as one goes through the reso- 
nance starting from the negative s-wave scattering length 
and passing to its positive values, the Gooper pairs of 
two atoms in two hyperfine states shrink in size and go 
over continuously into the Bose-Einstein condensate of 
diatomic molecules. In this respect, an important ques- 
tion arises as to what are the qualitative features allowing 
one to distinguish between the BGS and BEG states. 

In the low density limit, when the interaction is char- 
acterized by a single parameter, the s-wave scattering 
length a, the ^Sq pairing can be described on the base of 
the so-called regularized gap equation [2^ . In this model, 
the strength of the interaction is controlled only by the 
dimensionless parameter l/kpa. The range l/kpa <^ —1 
corresponds to the weak coupling BGS regime while the 
range l/kpa ^ 1 is related to the strong coupling BEG 
regime. The boundaries of the crossover region between 



the BGS and BEG states, according to Ref. [20|, are de- 
termined by the values l/kpa = ±1. The advantage 
of this regularized model is that all integrals in the gap 
equation and the particle number equation can be ex- 
pressed in terms of some special functions and, hence, 
the model is easily solvable. Then all quantities of in- 
terest can be calculated, in particular, at the boundaries 
l/kpa = ±1 of the crossover region and these values can 
be used as the reference values for the study of the BGS- 
BEG transition in a specific Fermi system at low density. 

Such a viewpoint was adopted in Refs. [H, for the 
study of a possible BGS-BEG crossover of neutron pairs 
in nuclear matter. There were considered a few prob- 
able characteristics of the crossover, such as the ratios 



ii-ms / dn, An/ spn [ij, llol , ^Ln/spn UM: ^nd the probabil- 
ity P{dn) for the pair neutrons to come close to each other 
within the relative distance dn [l3l. [isj (^rms, c^n, A„, £i?„ 
and /i„ being the r.m.s. radius of the Gooper pair, the 
mean average distance between neutrons, the neutron 
energy gap, the neutron Fermi kinetic energy, and the 
effective neutron chemical potential, respectively). It 
was shown that these criteria are closely related to each 
other and give approximately the same density range 
g-n/go ^ 10~^ — 10~^ for the domain, where the neu- 
tron pairs are strongly spatially correlated [l^l- Besides, 
in terms of these criteria the BGS-BEG crossover was 
found for a di-neutron condensate in symmetric nuclear 
matter with the density dependent contact interaction at 
neutron densities gn/ go ^ 10~^ [Tsj . 

However, it is necessary to note that in the case of nu- 
clear matter the assumption of the dilute gas limit may 
be justified only at very low densities gn/ go S 10~^ ]V^ . 
Hence, the use of the regularized model with the contact 
interaction for getting the reference values of various cri- 
teria can lead only to the qualitative results at the above 
densities. 

Thus, it is desirable to approach the problem from the 
different point of view and to look for some other possible 
characteristics for identifying the BGS and BEG states. 
In this respect, using the analogies with other low-density 
Fermi systems, such as, e.g., ultracold atomic gases, could 
be expedient. To that end, let us mention the idea to 
utilize the density-density correlations in the image of an 
expanding gas cloud to detect superfluid states in the 
system of fermionic atoms released from the trap [27| . 
This idea was explored in Ref. [2^ to study the alter- 
ations in the density-density correlations through the 
BGS-BEG crossover. It was learned that the density cor- 
relation function of two-component ultracold fermionic 
atomic gas with singlet pairing of fermions changes sign 
at low momentum transfer across the BGS-BEG tran- 
sition. This feature was considered as an unambiguous 
signature of the crossover. Here we would like to check 
this criterion with respect to the BGS-BEG crossover in 
a di-neutron condensate of nuclear matter. Besides, the 
calculations of Ref. [l^ , performed for the balanced pop- 
ulations of fermions of two species, will be extended to 
unbalanced populations of unlike fermions (for neutron- 
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proton systems, this means finite isospin asymmetry). In 
defining a primary characteristic for identifying tlic BCS- 
BEC transition, we adopt a viewpoint that tlie qualita- 
tive boundary between the BCS and BEC states occurs 
when the chemical potential reaches the zero value. Just 
this feature was considered as an indication of the BCS- 
BEC crossover in the earlier researches on the BCS-BEC 
transition in nuclear matter 0, Q . 



II. BASIC EQUATIONS 

The equation for the neutron energy gap in the 
pairing channel reads [1^ 



^"W = -^7Lnk,k)2^-^tanh^^. (1 



Here 



Enik) = VC2(k) + A2(k), $„(k) = e„(k) 



e„(k) and being the neutron single particle energy 
and neutron chemical potential, respectively. Further 
the single particle energy will be taken in the form cor- 
responding to the Skyrme effective interaction, e„ = 
fi^k'^ /2mn + Um with to„ and Un being the neutron effec- 
tive mass and neutron single particle potential, respec- 
tively. The quadratic on momentum term in the neutron 
single particle potential is already included into the ki- 
netic energy term and the remaining part Un is momen- 
tum independent and can be incorporated into the effec- 
tive neutron chemical potential /i„ = /^on ~ Un- Eq. ([T]) 
should be solved self-consistcntly with the equation for 
the neutron particle number density g„. 
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^ i?„(k)*^"^ 2T )-V ^ 



Let us introduce the neutron anomalous distribution 
function [2§| 

,,,, / + + \ A„(k) ^ Enjk) 

where (...) = Tig..., g being the density matrix of the 
system. Then, using Eq. one can represent Eq. ^ 
for the energy gap in the form 
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■V(k) + (1 - 2/k) ^(k, k')V'(k') = 2M„V(k). (3) 



In the limit of vanishing density, fk ^ 0, Eq. ([3|) goes 
over into the Schrodinger equation for a di-neutron bound 
state. The corresponding energy eigenvalue is equal to 
2/x„. The change of sign of the effective neutron chemical 
potential /x„ under decreasing density of nuclear matter 
signals the transition from the regime of large overlapping 



neutron Cooper pairs to the regime of non-overlapping 
di-neutron bound states. 

Let us consider the two-neutron density correlation 
function 

V{x, x') = (A7i(x)An(x')), A7i(x) = 7i(x) - n, (4) 



J(k'-k)x 4- 



Tkk' 



V 4^ 



Its general structure in the spatially uniform and 
isotropic case reads (s^l 



2?(x, x') = gnS{r) + gnD{r), r = x - x'. 



(5) 



The function D{r) is called the density correlation func- 
tion as well. We will just be interested in the behavior of 
the function D{r). Calculating the trace in Eq. (|4|) and 
going to the Fourier representation 



D{q) 



one can get 
where 



D{q)^I^,{q)^If{q), 



(6) 



T^-'Qn Jo 



■^"Qn Jo 



drr'^joirq) / dkk^ f{k)jo{rk) 

poo 

drr^jo{rq) I dk k'^%p{k)jo{rk) 



Here jo is the spherical Bessel function of the first kind 
of order zero. The functions If and represent the nor- 
mal and anomalous contributions to the neutron density 
correlation function. 

Analogously, we can consider the two-neutron spin cor- 
relation function 



>S^,(x,x') = (As^(x)As,(x')), As^(x) = s^(x) 



(7) 



E ^n,kcr('''/i)<TCT'an,ko-') 



2V 



where are the Pauli matrices. The general structure 
of the neutron spin correlation function for the spin un- 
polarized case is 

5^1. (x,x') = ^ 5^^ 6{r) + gn S^^{r). (8) 

Calculating the trace in Eq. (O, for the Fourier transform 
of the spin correlation function one gets 

S^^q) = S{q)S^,, S{q) = -\{If{q) + l^{q)). (9) 
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At zero momentum transfer, the normal If and anoma- 
lous I^, contributions read 



1 

T^^Qn Jo 
l4,iq = 0) = ^ 



dkk^fik), (10) 

dk k^1j?(k). 



At zero temperature, the zero-momentum transfer neu- 
tron spin correlation function satisfies the sum rule 

^(g = 0) = -i(//(0) + /v,(0))--J. (11) 



III. ^So NEUTRON ENERGY GAP AT ZERO 
TEMPERATURE 

Further we will take the pairing interaction in the ^5*0 
neutron pairing channel in the form of the effective den- 
sity dependent contact interaction 



y(ri, r2) = vog{g„, gp)S{ri - r2). 



(12) 



In the gap equation, the interaction ()12p should be sup- 
plemented with the cut-off energy Ec- The interaction 
strength vq is determined in such a way that the inter- 
action reproduces the low-energy neutron-neutron 
scattering phase shifts obtained with the Argonne vig po- 
tential. The isospin dependent factor g was constructed 
in Ref. [la| and reads 



9 = 91 +92, 



51 = 1 - (1 - a)7js — 



(13) 



arjn — 
\9o 



I3„ 



92 = m 



1 4- e 0.05 



1 + e 0.05 



Here a = (p„ — gp)/g is the isospin asymmetry parame- 
ter, kpn = (37r^gi„)^/'^ is the neutron Fermi momentum, 
Ps, l3m'ils','nniV2 sne the fitting parameters, adjusted to 
reproduce the pairing gap in symmetric nuclear matter 
and neutron matter obtained in the G-matrix calcula- 
tions with the Argonne vig potential with account of 
the medium polarization effects [l6j . The term gi con- 
tains the density dependence in a power law form as it 
was introduced in the contact interaction of Ref. [sij . 
Besides, the term gi contains the isospin dependence, 
built such as to reproduce the position and the magni- 
tude of the maximum of the energy gap in both sym- 
metric nuclear matter and neutron matter. The term 52 
is necessary to improve the results of the fitting proce- 
dure for the energy gap at the neutron Fermi momenta 
kpn > lfm~^. Concerning the term 52, let us make one 
more comment. In the original article p^ . this term is 
written as (72 = g2{kF), where the quantity A;f is not 
explicitly determined. There are two possibilities: (1) 
92 = 92{kF), where kp is parametrized in terms of the 



TABLE I: The parameters of the screened-II interaction [15| 
obtained from the fit to the neutron energy gap in symmetric 
nuclear matter and neutron matter, calculated for the Ar- 
gonne ?;i8 potential with account of the medium polarization 
effects. 



vo (MeVfm=*) 


(MeV) ri2 


■qs Ps Vn 


/3n 


-542 


40 0.8 


1.80 0.27 1.61 


0.122 



total density kp = (37r^p/2)^, and (2) 52 = 92{kpn) 
with kpn being the neutron Fermi momentum. There 
will be no difference between these two parametrizations 
for symmetric nuclear matter; however, the difference oc- 
curs for isospin asymmetric matter and will be most se- 
rious for neutron matter. The calculations of the energy 
gap in neutron matter show that the results of the fitting 
procedure, presented in Ref. [Tl|, correspond to the case 
92 = Q2(kpn). This point was confirmed in the commu- 
nication [32| . Therefore, the term 52 is, in fact, isovector 
and not isoscalar, as it is called in Ref. [lH|. Since the 
microscopically calculated neutron energy gap in sym- 
metric nuclear matter and neutron matter is well fitted 
with 92 = 92{kpn), one can use this form also for calcula- 
tions at intermediate asymmetries < a < 1. Besides, if 
to apply the isospin dependent contact interaction p^ . 
(|13p for the calculation of the pairing gaps in finite nu- 
clei [3^, one has to change a — *■ —a in the term gi and 
to use 32 = 92{kpp) {kpp being the proton Fermi momen- 
tum) in the proton ^Sq pairing channel. 

Let us now turn to the calculation of the neutron en- 
ergy gap in the ^Sq pairing channel that is the necessary 
preliminary step before the calculation of the correlation 
fimctions. The zero temperature gap equation with the 
effective contact interaction reads 



-vo9 



d^k 



/c_ <k<k4 



(27r)3 2Er,{k)' 



(14) 



In the gap equation, the cut-off momenta k± are de- 
fined from the requirement that the quasiparticle energy 
should be less than the cut-off energy Ec, En{k) < Ec- If 
k- becomes imaginary, one should set fc_ = 0. Just this 
choice of the cut-off parameters is often used in the HFB 
calculations. Eq. (fTi|) should be solved self-consistently 
with Eq. (12) for the neutron particle number density. In 
numerical calculations, we use the screened-II interaction 
of Ref. [1^ with the parameters given in Table 1. The 
single particle energy £n(fc) is taken for the SLy4 Skyrme 
effective interaction. 

In Fig. [1^, the neutron energy gap is plotted as a func- 
tion of the neutron Fermi momentum at zero tempera- 
ture and different isospin asymmetries. With increasing 
asymmetry, the magnitude of the energy gap decreases 
and the position of its maximum is shifted to the larger 
neutron Fermi momentum. In Fig. [TJd, the ratio of the 
effective neutron chemical potential ^„ to the neutron 
Fermi kinetic energy epyi = 7i^A;|i„/2m„ is depicted as 
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FIG. 1: (Color online) (a) The zero temperature neutron en- 
ergy gap and (b) the ratio /i„/e„ of the effective neutron 
chemical potential to the neutron Fermi kinetic energy as 
functions of the neutron Fermi momentum at different isospin 
asymmetries for the screened-II contact interaction. 

a function of the neutron Fermi momentum. This ratio 
is indicative of the importance of the di-neutron corre- 
lations, since in the absence of pairing correlations it is 
equal to unity, otherwise it is less than unity. It is seen 
that the di-neutron correlations are most strong in sym- 
metric nuclear matter and with increasing isospin asym- 
metry their strength is diminished. In accordance with 
Eq. we wiU identify the BCS and BEC states of 
the pair condensate by the positive (BCS regime) and 
negative (BEC regime) values of the effective chemi- 
cal potential /i„. For symmetric nuclear matter, a di- 
neutron BEC state is realized at the finite density range 
0.09 fm^^ < kpn < 0.25 fm^^ and surrounded on both 
sides by the BCS state. For neutron matter, the self- 
consistent equations have no at all solutions with the 
negative chemical potential /i„ and the BEC state is not 
formed. For asymmetric nuclear matter, the BEC state is 
quite sensitive to the isospin asymmetry and is quenched 
at the critical value ac ~ 0.26. Note that, as was clarified 



FIG. 2: (Color online) The normal f{k) and anomalous ^p{k) 
distribution functions of neutrons in symmetric nuclear mat- 
ter at different neutron Fermi momenta. 



in Ref. [l6j . the medium polarization acts differently on 
the di-neutron correlations at low-density region, increas- 
ing the neutron energy gap in symmetric nuclear matter 
and decreasing it in neutron matter. 

Fig. [2] shows the evolution of the neutron normal f{k) 
and anomalous Tjj{k) distribution functions with density 
in symmetric nuclear matter. In the weak coupling BCS 
regime (kpn = 0.8 fm"^), the normal distribution func- 
tion resembles the Fermi step distribution function with 
the edge near the neutron Fermi momentum. With 
decreasing density, the step-like structure is gradually 
washed out and in the region of the strong coupling BEC 
regime (kpn = 0.2 fm~^) the normal distribution func- 
tion smoothly changes with the momentum. However, 
under further decreasing density, we go back to the BCS 
regime {kpn = 0.02 fm~^) and the step-like structure be- 
comes more and more pronounced. The anomalous dis- 
tribution function ?/'(^) in the BCS regime has a well de- 
veloped peak, indicating that pairing exists only near the 
Fermi surface. In the strong coupling BEC regime ip{k) 
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FIG. 3: (Color online) (a) Density D{q) and (b) spin S{q) 
correlation functions as functions of the momentum trans- 
fer at different neutron Fermi momenta in symmetric nuclear 
matter. 



FIG. 4: (Color online) Zero-momentum transfer density cor- 
relation function D{q = 0) and its normal //(g — 0) and 
anomalous I^{q = 0) contributions as functions of the neu- 
tron Fermi momentum in symmetric nuclear matter. 



exhibits a smooth variation typical of the wave function 
of a bound state. At ultra-low densities, the spike in 
the momentum dependence of the anomalous distribu- 
tion function appears again. 

IV. DENSITY AND SPIN CORRELATION 

FUNCTIONS FOR A DI-NEUTRON 
CONDENSATE AT ZERO TEMPERATURE 

A. Symmetric nuclear matter 

After the neutron energy gap and the effective neutron 
chemical potential have been found, the density D(q) 
and spin S{q) correlation functions can be obtained di- 
rectly from Eqs. ^ and Fig. [3^ shows the den- 
sity correlation function as a function of the momen- 
tum transfer for a set of the neutron Fermi momenta 
in symmetric nuclear matter. Let us begin with the 
consideration of the BCS regime at the moderate den- 



sities {kpn = 0.8 fm~ ), where at low momentum trans- 
fer D(q) is negative. Under decreasing density {kpn = 
0.6, 0.4fm~^), the zero-momentum transfer density cor- 
relation function D{q = 0) at certain neutron Fermi mo- 
mentum changes sign from negative to positive. In the 
BEC regime (kpn = 0.2 fm~^), the correlation function 
D{q) is everywhere positive. If to decrease density fur- 
ther (kpn = 0.07, 0.02 fm~^), the zero-momentum trans- 
fer value D(q = 0) changes sign again, now from posi- 
tive to negative, and in the ultra-low density BCS regime 
the function D(q) has the negative sign at low momen- 
tum transfer as it was at the moderate densities. The 
change of sign of the density correlation function at low 
momentum transfer across the BCS-BEC transition re- 
gion was considered in Ref. [1^ as a signature of the 
crossover. We sec that this feature is qualitatively well 
reproduced for a di-ncutron condensate in symmetric nu- 
clear matter. At the same time, the spin correlation func- 
tion, as can be seen from Fig. fluently evolves across 
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FIG. 5: (Color online) Density correlation function D{q) as 
a function of the momentum at different neutron Fermi mo- 
menta and isospin asymmetries. 



the BCS-BEC transition region without any qualitative 
change. In order to establish the boundaries of the den- 
sity interval for the BEC state on the basis of the for- 
mulated criterion, let us consider the density correlation 
function at zero momentum transfer, depicted in Fig. |4l 
Under decreasing density, the density correlation func- 
tion D{q = 0) twice changes sign, first, at the point of 
the BCS-to-BEC transition and then at the point of the 
reverse BEC-to-BCS crossover, that corresponds to the 
earlier consideration. The BEC state is realized in the 
range 0.03 fm~^ < kpn < 0.5 fm~^ and this interval is 
wider than that determined on the basis of the negative 
values of the effective neutron chemical potential. 

The normal If{q = 0) and anomalous I^{q — 0) contri- 
butions to D{q ~ 0), given by Eq. p0|) . are shown in the 
bottom panel of Fig. (4] This figure qualitatively explains 
why the density correlation function D{q = 0) changes 
sign at the boundary between the BCS and BEC states. 
In the weak couphng BCS regime {kpn = 0.8 fm~^), the 
anomalous distribution function ip^k) is sharply peaked 
near k = kpn, as shown in Fig[2l and, hence, the anoma- 
lous contribution to D{q = 0) is much smaller than the 
normal one. Therefore, D{q = 0) < on the BCS side 
from the transition boundary. However, in the strong 
coupling BEC regime {kpn = 0.2 fm~^), the normal dis- 
tribution function f{k) is to a considerable extent de- 
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FIG. 6: (Color online) Zero-momentum transfer density cor- 
relation function D(g = 0) as a function of the neutron Fermi 
momentum at different isospin asymmetries. 



pletcd and the anomalous contribution dominates over 
the normal one. Hence, D((i = 0) > on the BEC side 
from the transition boundary, and at some critical den- 
sity between the BCS and BEC regimes the density cor- 
relation function £'((7 = 0) should change the sign. Note 
also that, according to the sum rule in Eq. pT|) . the nor- 
mal and anomalous contributions satisfy the relationship 
//(0) + /^(0) = l. 

As seen from the figure, the contribution of the anoma- 
lous part is the largest at the neutron Fermi momen- 
tum kpn ^ 0.2 fm~^. At the corresponding densities 
Qn/ Qq ^ 10^'^ the BEC di-neutron state can be consid- 
ered as well formed. Just at these densities the BEC of 
bound neutron pairs in symmetric nuclear matter was 
found in Ref. [l^ with the help of other quantitative 
characteristics, mentioned in Introduction. 



B. Asymmetric nuclear matter 

Now wc will present the results of the numerical de- 
termination of the correlation functions at finite isospin 
asymmetry. The evolution of the density correlation 
function D[q) with the isospin asymmetry is shown in 
Fig. m If to identify the BCS-BEC transition by the 
change of sign of the density correlation function at low 
momentum transfer, it is seen that with increasing asym- 
metry the density interval, where the BEC state exists, 
is contracted. At strong enough isospin asymmetry only 
the BCS state is realized, as can be seen from the limiting 
case a = 1, corresponding to neutron matter. 

In order to find an estimate for the critical isospin 
asymmetry, at which the BEC state disappears, let us 
consider the zero-momentum transfer density correlation 
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FIG. 7: (Color online) Spin correlation function S{q) as a 
function of the momentum at different neutron Fermi mo- 
menta and isospin asymmetries. 

function D{q = 0) as a function of the neutron Fermi 
momentum at different asymmetries (Fig. [5]). For large 
isospin asymmetry a > af ^ 0.9 the function D{q = 0) 
is always negative that corresponds to the BCS state for 
all densities where a di-neutron condensate exists. This 
critical asymmetry is much larger than ac ~ 0.26, found 
earlier on the basis of the change of sign of the chemical 
potential fin- From that one can infer that the crite- 
rion of the BCS-BEC crossover, based on the change of 
sign of the density correlation function, can be trusted 
only at small isospin asymmetry. This conclusion qual- 
itatively agrees with that of Ref. [U, where the den- 
sity and spin correlation functions were calculated for a 



neutron-proton condensate in dilute asymmetric nuclear 
matter. 

Fig. [7] shows the evolution of the spin correlation func- 
tion S{q) with the isospin asymmetry. It is seen that the 
zero-momentum transfer value S{q = 0) satisfies the sum 
rule (jlip . independently of the isospin asymmetry. The 
spin correlation function S{q) smoothly varies between 
its limiting values without any qualitative change. 

In summary, the density and spin correlation func- 
tions have been calculated for a di-neutron condensate 
in asymmetric nuclear matter with the aim of finding 
the possible signatures of the BCS-BEC crossover. As 
the primary characteristic for identifying the transition, 
it is accepted that the qualitative boundary between the 
BCS and BEC states occurs when the effective neutron 
chemical potential reaches the zero value. It has been 
shown that the zero-momentum transfer spin correlation 
fimction satisfies the sum rule at zero temperature. In 
symmetric nuclear matter, the density correlation func- 
tion changes sign at low momentum transfer across the 
BCS-BEC transition and this feature can be considered 
as a signature of the crossover. This result qualitatively 
agrees with that obtained for two-component ultracold 
fermionic atomic gas with equal densities of fermions of 
different species The self-consistent calculations of 
the ^5*0 energy gap in a di-neutron condensate and the ef- 
fective neutron chemical potential /x„ show that the BEC 
state (/x„ < 0) is quite sensitive to isospin asymmetry 
and does not appear at a > etc ~ 0.26. This is a reflec- 
tion of the fact that the medium polarization effects act 
differently on di-neutron correlations at small and strong 
isospin asymmetry, enhancing them in symmetric nuclear 
matter and suppressing them in neutron matter. At finite 
isospin asymmetry, the criterion, based on the change of 
sign of the density correlation function, gives too large 
value for the critical asymmetry aj? ^ 0.9, at which the 
BEC state is quenched. Therefore, it can be trusted for 
the description of the density-driven BCS-BEC crossover 
of neutron pairs only at small isospin asymmetry. This 
result generalizes the conclusion of Ref. [28| , in which 
the change of sign of the density correlation function 
at low momentum transfer in two-component quantum 
fermionic atomic gas with the balanced populations of 
fermions of different species was considered as an unam- 
biguous signature of the BCS-BEC transition. 
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